This project seeks to estimate sport fish harvest and releases of rockfish in Alaska waters by improving on the Howard et al. (2020) methods and expand the time series back to 1977 when the statewide harvest survey (SWHS) was first implemented. This is essentially a Bayesian version of the Howard methods that allows for more appropriate and defensible sharing of information between areas, handles missing data in a more appropriate manor, accurately propagates uncertainty throughout the estimation procedure and replaces the Howard decision tree approach to low sample sizes with a hierarchical model. The methods and results for generating harvest estimates are generally consistent between the Bayesian model and the Howard methods. Harvest estimates are consistent with Howard estimates during contemporary times, but may differ based on more appropriate weighting of SWHS and logbook data, including estimating and correcting bias in the SWHS data.

Table 1. Summary of key improvements in reconstructiing sport fish removals of rockfish using the Bayesian model as compared to the Howard et al. (2020) methods.
Issue Howard Bayes
Time series 1999 - present 1977 - present
Bias in SWHS Not explicitly dealt with. Relies on logbook data and ratios of guided/unguided from SWHS data to estimate unguided releases and harvests. Explicitly estimates bias in SWHS harvest and release estimates based on logbook data.
Species composition of releases Assumes that species composition of releases is equal to that of the harvest, which is not evident in the logbook data. Recognizes different release probabilities by species / species assemblage and estimates it from logbook data and bias corrected SWHS data
Sample size limitations Uses sample size threshholds such that when areas fall below those threshholds values are borrowed from nearby areas. Uses a hierarchichacal modelling approach that shares information between areas in the same region. Thus all data is used, even with small sample sizes. This is a more sound method that avoids assumptions and uses all of the data.
Error propogation Error is propogated when variance estimates are available, but there is uncertainty associated with borrowing values from nearby areas, or the assumption of species compositions being identical in harvest and releases, are not dealt with. By breaking the assumption that species composition is equal between harvests and releases, uncertainty in the release estimates is more reflective of the fishery. Furthermore, the hyerarchichal approach more accurately captures uncertainy within and between areas within a region.
**Figure 1.**- Data sources for estimating rockfish harvests and releases in ADF&G commercial fisheries management units. Note that initial rockfish harvest estimates were not differentiated into species assemblage or species until 1998 when logbooks began differentiating by pelagic and non-pelagic. Logbooks began to collect data on yelloweye beginning in 2006. Port sampling programs to gather data on species composition of harvests began in 1996 in Southcentral and Kodiak and in 2006 in Southeast.

Figure 1.- Data sources for estimating rockfish harvests and releases in ADF&G commercial fisheries management units. Note that initial rockfish harvest estimates were not differentiated into species assemblage or species until 1998 when logbooks began differentiating by pelagic and non-pelagic. Logbooks began to collect data on yelloweye beginning in 2006. Port sampling programs to gather data on species composition of harvests began in 1996 in Southcentral and Kodiak and in 2006 in Southeast.


**Figure 2.**- SWHS harvest (left) and release (right) estimates from guided trips (x-axis) versus repoted harvests from charter logbooks (y-axis).

Figure 2.- SWHS harvest (left) and release (right) estimates from guided trips (x-axis) versus repoted harvests from charter logbooks (y-axis).


Results

**Figure X.**- Rhat values and proportion of parameters that converged (Rhat < 1.1.)

Figure X.- Rhat values and proportion of parameters that converged (Rhat < 1.1.)

Estimate comparison

Since previous estimates of rockfish harvest have been produced these first 3 graphs will be used to show how the modeled estimates compare to the estimates produced earlier. For total rockfish the estimates are in general agreement although differences are noted. These estimates should be more reliable because they include both SWHS and guide logbook data, handle variance more appropriately, use hierarchical distributions when data is missing, directly consider observation error and are produced using reproducible research.

**Figure 3.**- Total rockfish harvests 1996-2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 3.- Total rockfish harvests 1996-2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



**Figure 3.**- Total rockfish releases 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 3.- Total rockfish releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.


Notes from Adam: When looking at only black rockfish the most significant differences are for the Prince William Sound Inside area. I did not spend a great deal of time tracking this down although it looks like the previous version used bad values for \(P_{(black)ayu}\) for at least unguided anglers. For the moment I would ignore the results for BSIA and SOKO2SAP. I think it is possible to give approximate values for these areas but it will require a little more coding which I have yet to do.

**Figure 4.**- Black rockfish harvests 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 4.- Black rockfish harvests 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.


And black rockfish releases…

**Figure 5.**- Black rockfish releases 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 5.- Black rockfish releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.





**Figure 6.**- Yellow rockfish harvests 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 6.- Yellow rockfish harvests 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



**Figure 7.**- Yellow rockfish releases 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 7.- Yellow rockfish releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.





**Figure 8.**- DSR rockfish (excluding yelloweye) harvests 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 8.- DSR rockfish (excluding yelloweye) harvests 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



**Figure 9.**- DSR rockfish releases (including yelloweye) 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 9.- DSR rockfish releases (including yelloweye) 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



**Figure 11.**- Slope rockfish harvests 1996-2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 11.- Slope rockfish harvests 1996-2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



**Figure 12.**- Slope rockfish releases 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 12.- Slope rockfish releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

**Figure 13.**- *Other* rockfish (non-pelagic and non-yelloweye) harvests 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 13.- Other rockfish (non-pelagic and non-yelloweye) harvests 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



**Figure 14.**- *Other* rockfish (non-pelagic and non-yelloweye) releases 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 14.- Other rockfish (non-pelagic and non-yelloweye) releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



Total Biomass Removal Estimates

**Figure 15.**- Black rockfish estimated total removals in lbs in 1996--2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Figure 15.- Black rockfish estimated total removals in lbs in 1996–2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.



**Figure 16.**- Yellow rockfish estimated total removals in lbs in 1996--2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Figure 16.- Yellow rockfish estimated total removals in lbs in 1996–2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

**Figure 17.**- Pelagic rockfish estimated total removals in lbs in 1996--2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Figure 17.- Pelagic rockfish estimated total removals in lbs in 1996–2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.


**Figure 18.**- Non-yelloweye, demersal shelf rockfish estimated total removals in lbs in 1996-2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Figure 18.- Non-yelloweye, demersal shelf rockfish estimated total removals in lbs in 1996-2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.


**Figure 19.**- Slope rockfish estimated total removals in lbs in 1996-2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Figure 19.- Slope rockfish estimated total removals in lbs in 1996-2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.


Model fit

Logbook residuals

**Figure 20.**- Residuals from logbook harvests.

Figure 20.- Residuals from logbook harvests.


SWHS residuals

**Figure 21.**- Residuals from SWHS harvests.

Figure 21.- Residuals from SWHS harvests.



**Figure 22.**- Residual of SWHS releases.

Figure 22.- Residual of SWHS releases.

Parameter estimates

P(Charter)

These histograms show the posterior distribution of the mean percent of rockfish harvested by the charter fleet.

**Figure 23.**- Mean percent of harvest by charter anglers.

Figure 23.- Mean percent of harvest by charter anglers.


When considered annually we see the percent of rockfish harvested by the charter fleet follows our data fairly well although the model smooths out the changes and we just do not have much information about this ratio. Prior to 2011 the percent charter is confounded with SWHS bias and should be mostly discounted.

**Figure 24.**- Annual estimates of the percent of harvest by charter anglers for 16 commerical fishing manamgent areas, 1996-2023.

Figure 24.- Annual estimates of the percent of harvest by charter anglers for 16 commerical fishing manamgent areas, 1996-2023.

P(Harvest)

These plots show the fitted logistic line to the proportion of caught rockfish that are harvested. These estimates are used for hindcasting catch estimates based on the harvest data in early years when catch estimates are unavailable.


**Figure 25.**- Annual proportion of pelagic rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available.

Figure 25.- Annual proportion of pelagic rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available.


**Figure 26.**- Annual proportion of yelloweye rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available.

Figure 26.- Annual proportion of yelloweye rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available.


**Figure 27.**- Annual proportion of non-pelagic, non-yelloweye rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note, that this is not estimated for Southeast areas because non=pelagics are divided between DSR (including yelloweye) and Slope species.

Figure 27.- Annual proportion of non-pelagic, non-yelloweye rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note, that this is not estimated for Southeast areas because non=pelagics are divided between DSR (including yelloweye) and Slope species.


**Figure 28.**- Annual proportion of DSR rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note that the observed logbook data is for all non-pelagic, non-yelloweye fish.

Figure 28.- Annual proportion of DSR rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note that the observed logbook data is for all non-pelagic, non-yelloweye fish.


**Figure 29.**- Annual proportion of Slope rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note that the observed logbook data is for all non-pelagic, non-yelloweye fish.

Figure 29.- Annual proportion of Slope rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note that the observed logbook data is for all non-pelagic, non-yelloweye fish.

SWHS bias

Figure 23 shows the mean estimate for SWHS bias in harvests and releases. Cook Inlet, North Gulf Coast and North Southeast Inside all look pretty good while most other areas have substantial bias. Prince William Sound Inside has the largest bias. Bias in release estimates is substantial and whereas the SWHS appears to underestimate harvests, it appears to greatly overestimates releases by a factor of 2 or more in most areas as derived from logbook reported releases.

**Figure 30.**- Mean SWHS bias for harvests and catches. Note that a bias < 1 indicates that the SWHS *underestimates* the true value and bias > 1 indicates the survey *overestimates* the true value.

Figure 30.- Mean SWHS bias for harvests and catches. Note that a bias < 1 indicates that the SWHS underestimates the true value and bias > 1 indicates the survey overestimates the true value.


Our estimates of SWHS harvest bias track observations fairly well when he have guided harvest estimates. The estimates of release bias in the SWHS data track observed patterns to an extent, but appear to smooth these more volatile disagreements with the logbook data. Adam postulated in his initial start on this that some of this could be the result of the estimates of the proportion guided. This value was not modelled with a trend and thus applies a constant estimate when hindcasting. Data on these relationships could greatly improve this model.

**Figure 31.**- Annual estimates of SWHS bias in harvests and releases for 16 commerical fishing manamgent areas, 1996-2023. Note that a bias < 1 indicates that the SWHS *underestimates* the true value and bias > 1 indicates the survey *overestimates* the true value.

Figure 31.- Annual estimates of SWHS bias in harvests and releases for 16 commerical fishing manamgent areas, 1996-2023. Note that a bias < 1 indicates that the SWHS underestimates the true value and bias > 1 indicates the survey overestimates the true value.

P(pelagic)

We model the percentage of pelagic rockfish in the harvest because we have the information for charter anglers (via logbooks) starting in 1998. Other than looking at the model estimates you can use Figure 25 to compare the two data streams for pelagic rockfish harvest. In general they are in agreement with major exceptions in Price William Sound inside, Prince William Sound outside (early in the time series) and South Southeast inside.

**Figure 32.**- Annual estimates of the percent of the sport harvest that was pelagic rockfish for 16 commerical fishing manamgent areas, 1996-2023.

Figure 32.- Annual estimates of the percent of the sport harvest that was pelagic rockfish for 16 commerical fishing manamgent areas, 1996-2023.

P(black|pelagic)

Note that in Southeast Alaska we only have composition data starting in 2006. Tania dug up old SE data, but it did not provide any useful data for species apportionment. For the most part, P(black|pelagic) is relatively constant across areas, with the exception of Cook Inlet and NSEI in Southeast AK. It may be worth discussing whether the shifts in those areas is a result of improved or changing species identification rather than actual shift in the species composition of the catch.

**Figure 33.**- Annual estimates of the percent of the sport harvest of pelagic rockfish that were black rockfish for 16 commerical fishing manamgent areas, 1996-2023. Kodiak panels include data from a hydroacoustic survey and the proportion of pelagic rockfish that are black in those areas (red) and the adjusted proportions based on obseved harvests for charter (blue) and private (cyan) users.

Figure 33.- Annual estimates of the percent of the sport harvest of pelagic rockfish that were black rockfish for 16 commerical fishing manamgent areas, 1996-2023. Kodiak panels include data from a hydroacoustic survey and the proportion of pelagic rockfish that are black in those areas (red) and the adjusted proportions based on obseved harvests for charter (blue) and private (cyan) users.

P(yelloweye|non-pelagic / yelloweye|DSR)

**Figure 34.**- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were yelloweye rockfish for 16 commerical fishing manamgent areas, 1996-2023. Note that P(yelloweye) is the the proportion relative to non-pelagics for Central and Kodiak areas but is relative to DSR for Southeast areas.

Figure 34.- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were yelloweye rockfish for 16 commerical fishing manamgent areas, 1996-2023. Note that P(yelloweye) is the the proportion relative to non-pelagics for Central and Kodiak areas but is relative to DSR for Southeast areas.

P(DSR|non-pelagic)

**Figure 35.**- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were DSR rockfish for 6 Southeast commerical fishing manamgent areas, 1996-2023.

Figure 35.- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were DSR rockfish for 6 Southeast commerical fishing manamgent areas, 1996-2023.

P(slope|non-pelagic)

**Figure 36.**- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were slope rockfish for 6 southeast commerical fishing manamgent areas, 1996-2023. Note that P(yelloweye) is the the proportion relative to non-pelagics for Central and Kodiak areas but is relative to DSR for Southeast areas.

Figure 36.- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were slope rockfish for 6 southeast commerical fishing manamgent areas, 1996-2023. Note that P(yelloweye) is the the proportion relative to non-pelagics for Central and Kodiak areas but is relative to DSR for Southeast areas.



P(slope|non-pelagic & non-yellowye) For release estimates

## NULL



Weight Fits

**Figure 37.**- Mean weights of black rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Figure 37.- Mean weights of black rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.


**Figure 38.**- Mean weights of yelloweye rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Figure 38.- Mean weights of yelloweye rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.


**Figure 39.**- Mean weights of non-black, pelagic rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Figure 39.- Mean weights of non-black, pelagic rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.


**Figure 40.**- Mean weights of non-yelloweye, demersal shelf rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Figure 40.- Mean weights of non-yelloweye, demersal shelf rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.


**Figure 41.**- Mean weights of slope rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Figure 41.- Mean weights of slope rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.


Summary of unconverged parameters:

The proportion of unconverged species composition proportional estimates that are near the bounds of 0 or 1 is 0.6879433

The proportion of unconverged true parameters (defined as Rhat > 1.01) that are also in excess of an Rhat of 1.1 is 0.09 while 0.91 are less than 1.1.

Table Xa. Species proportion derived quantities that are not with 0.05 of the 0 and 1 bounds and are not converged relative to an Rhat of 1.01.

variable

area_n

year

user

mean

median

lower

upper

Rhat

sd

skew

area

p_pelagic

13

1

1

0.228

0.205

0.053

0.514

1.071

0.123

0.187

NSEI

p_pelagic

13

1

2

0.090

0.074

0.017

0.253

1.052

0.062

0.258

NSEI

p_pelagic

13

10

1

0.258

0.240

0.079

0.535

1.038

0.121

0.149

NSEI

p_pelagic

13

10

2

0.102

0.089

0.025

0.260

1.064

0.061

0.212

NSEI

p_pelagic

13

11

1

0.267

0.249

0.084

0.551

1.025

0.121

0.149

NSEI

p_pelagic

13

11

2

0.107

0.093

0.026

0.273

1.027

0.065

0.217

NSEI

p_pelagic

13

12

1

0.274

0.259

0.090

0.545

1.032

0.119

0.126

NSEI

p_pelagic

13

12

2

0.108

0.095

0.029

0.267

1.022

0.062

0.211

NSEI

p_pelagic

13

13

1

0.280

0.265

0.092

0.551

1.025

0.120

0.125

NSEI

p_pelagic

13

13

2

0.112

0.099

0.030

0.267

1.018

0.064

0.203

NSEI

p_pelagic

13

14

1

0.275

0.262

0.096

0.520

1.024

0.109

0.119

NSEI

p_pelagic

13

14

2

0.114

0.101

0.032

0.273

1.017

0.062

0.209

NSEI

p_pelagic

13

15

1

0.250

0.235

0.089

0.488

1.014

0.102

0.147

NSEI

p_pelagic

13

15

2

0.110

0.099

0.034

0.247

1.021

0.056

0.196

NSEI

p_pelagic

13

16

1

0.284

0.269

0.110

0.522

1.011

0.110

0.136

NSEI

p_pelagic

13

16

2

0.123

0.109

0.034

0.287

1.018

0.065

0.214

NSEI

p_pelagic

13

17

1

0.303

0.293

0.122

0.551

1.016

0.110

0.091

NSEI

p_pelagic

13

17

2

0.130

0.119

0.040

0.288

1.016

0.064

0.171

NSEI

p_pelagic

13

2

1

0.234

0.214

0.058

0.536

1.064

0.125

0.160

NSEI

p_pelagic

13

2

2

0.091

0.074

0.018

0.245

1.070

0.063

0.271

NSEI

p_pelagic

13

3

1

0.235

0.215

0.057

0.516

1.060

0.121

0.165

NSEI

p_pelagic

13

3

2

0.091

0.077

0.017

0.249

1.047

0.061

0.228

NSEI

p_pelagic

13

4

1

0.238

0.217

0.056

0.528

1.050

0.124

0.169

NSEI

p_pelagic

13

4

2

0.092

0.078

0.018

0.255

1.052

0.061

0.230

NSEI

p_pelagic

13

5

1

0.240

0.218

0.061

0.521

1.054

0.122

0.180

NSEI

p_pelagic

13

5

2

0.093

0.078

0.021

0.240

1.034

0.060

0.251

NSEI

p_pelagic

13

6

1

0.240

0.220

0.065

0.526

1.051

0.121

0.166

NSEI

p_pelagic

13

6

2

0.095

0.080

0.022

0.250

1.052

0.060

0.249

NSEI

p_pelagic

13

7

1

0.245

0.226

0.069

0.529

1.039

0.120

0.158

NSEI

p_pelagic

13

7

2

0.098

0.084

0.020

0.262

1.035

0.062

0.225

NSEI

p_pelagic

13

8

1

0.250

0.229

0.074

0.521

1.030

0.119

0.176

NSEI

p_pelagic

13

8

2

0.098

0.084

0.023

0.256

1.037

0.062

0.228

NSEI

p_pelagic

13

9

1

0.257

0.237

0.076

0.534

1.039

0.121

0.166

NSEI

p_pelagic

13

9

2

0.101

0.087

0.024

0.260

1.028

0.063

0.222

NSEI

p_pelagic

16

1

2

0.373

0.365

0.133

0.653

1.012

0.137

0.058

SSEO

p_pelagic

16

10

1

0.399

0.394

0.157

0.671

1.012

0.135

0.037

SSEO

p_pelagic

16

2

2

0.370

0.361

0.135

0.649

1.011

0.135

0.067

SSEO

p_pelagic

16

4

1

0.395

0.387

0.142

0.680

1.017

0.140

0.057

SSEO

p_pelagic

16

5

1

0.392

0.384

0.149

0.671

1.012

0.138

0.058

SSEO

p_pelagic

16

7

2

0.375

0.366

0.145

0.647

1.015

0.134

0.067

SSEO

p_pelagic

3

10

1

0.557

0.562

0.296

0.796

1.012

0.131

-0.038

PWSI

p_pelagic

3

7

1

0.557

0.565

0.279

0.790

1.011

0.132

-0.061

PWSI

p_pelagic

3

9

1

0.555

0.561

0.301

0.791

1.014

0.129

-0.047

PWSI

p_black

8

40

1

0.480

0.415

0.342

0.963

1.017

0.184

0.352

afognak

Table Xb. True parameters with an Rhat in excess of 1.1

variable

area_n

year/species

mean

median

lower

upper

Rhat

sd

skew

year

area

beta2_pH

11

4

-4.357

-3.381

-14.546

-0.638

1.130

3.591

-0.272

CSEO

beta2_pH

12

4

-4.917

-3.916

-14.520

-1.023

1.131

3.617

-0.277

EWYKT

beta2_pH

13

4

-4.507

-3.455

-14.779

-0.779

1.148

3.630

-0.290

NSEI

beta0_pelagic

13

-1.616

-1.498

-3.306

-0.381

1.193

0.780

-0.151

NSEI

beta1_pelagic

13

3.981

3.588

1.975

7.127

1.245

1.479

0.266

NSEI

beta2_pelagic

13

0.133

0.123

0.056

0.281

1.152

0.061

0.165

NSEI

beta2_pH

14

4

-4.556

-3.533

-14.529

-0.971

1.154

3.488

-0.293

NSEO

beta2_pH

15

4

-4.187

-3.299

-14.037

-0.754

1.318

3.357

-0.264

SSEI

beta2_pH

16

4

-4.655

-3.633

-14.399

-1.281

1.146

3.353

-0.305

SSEO

mu_beta2_pH

3

4

-3.871

-3.038

-13.508

-0.307

1.217

3.281

-0.254

PWSI

beta4_pelagic

5

0.941

0.592

-6.374

12.051

1.120

7.695

0.045

BSAI

lambda_H

15

0.021

0.014

0.003

0.065

1.274

0.099

0.071

SSEI

Table Xc. Summary statistics describing the characteristics of the derived harvest and release estimates that are not converged relative to an Rhat of 1.01.

Statistic

non.BSAI.SOKO.Harv

BSAI.SOKO.Harv

non.BSAI.SOKO.Rel

BSAI.SOKO.Rel

n unconverged at Rhat = 1.01

267.00

189.00

873.00

309.00

Prop. of total

0.59

0.41

0.74

0.26

Prop. median less than 1

0.30

0.24

0.33

0.48

Prop. median less than 10

0.46

0.63

0.54

0.81

Prop. median less than 50

0.58

0.84

0.68

0.87

Median of the median values

22.05

5.50

6.00

1.32

Mean of the median values

280.72

48.72

230.87

63.74

Mean skew of posteriors

0.26

0.30

0.27

0.24

Table Xd. Summary statistics describing the characteristics of the derived harvest and release estimates that are not converged relative to an Rhat of 1.1.

Statistic

non.BSAI.SOKO.Harv

BSAI.SOKO.Harv

non.BSAI.SOKO.Rel

BSAI.SOKO.Rel

n unconverged at Rhat = 1.1

6.00

6.00

32.00

41.00

Prop. of total

0.50

0.50

0.44

0.56

Prop. median less than 1

0.50

1.00

0.69

0.63

Prop. median less than 10

0.67

1.00

0.75

0.95

Prop. median less than 50

1.00

1.00

0.91

0.95

Median of the median values

3.95

0.15

0.01

0.09

Mean of the median values

10.47

0.15

18.58

13.48

Mean skew of posteriors

0.19

0.20

0.13

0.17

Table Xe. All parameters and derived quantities that are not converged relative to an Rhat of 1.1.

variable

Freq

beta0_pelagic

1

beta1_pelagic

1

beta2_pelagic

1

beta2_pH

6

beta4_pelagic

1

Bp_ayg

1

Hb_ay

1

Hb_ayg

1

Ho_ay

2

Ho_ayg

3

Ho_ayu

4

Hp_ay

1

Hp_ayg

1

Hs_ayu

1

Hy_ayg

2

lambda_H

1

mu_beta2_pH

1

R_ayg

3

R_ayu

5

Rb_ayg

3

Rb_ayu

1

Rd_ay

2

Rd_ayg

3

Rd_ayu

1

Ro_ay

7

Ro_ayg

10

Ro_ayu

17

Rp_ayg

3

Rp_ayu

1

Rs_ay

1

Rs_ayg

2

Rs_ayu

1

Ry_ay

14

Ry_ayg

4

Ry_ayu

20

Parameter estimates:

Summary Table of Parameter Estimates
Parameter mean sd Lower_CI Median Upper_CI
mu_bc_H[1] -0.019 0.069 -0.151 -0.020 0.122
mu_bc_H[2] 0.075 0.053 -0.021 0.072 0.184
mu_bc_H[3] -0.270 0.068 -0.397 -0.272 -0.131
mu_bc_H[4] -0.923 0.181 -1.279 -0.917 -0.577
mu_bc_H[5] 1.385 1.203 0.120 1.106 4.447
mu_bc_H[6] -2.402 0.306 -2.992 -2.408 -1.808
mu_bc_H[7] -0.452 0.104 -0.662 -0.449 -0.255
mu_bc_H[8] 0.247 0.339 -0.351 0.223 0.983
mu_bc_H[9] -0.248 0.126 -0.494 -0.254 0.011
mu_bc_H[10] 0.059 0.066 -0.062 0.057 0.197
mu_bc_H[11] -0.123 0.037 -0.198 -0.123 -0.051
mu_bc_H[12] -0.255 0.105 -0.478 -0.250 -0.056
mu_bc_H[13] -0.158 0.071 -0.295 -0.160 -0.015
mu_bc_H[14] -0.305 0.089 -0.493 -0.302 -0.137
mu_bc_H[15] -0.334 0.047 -0.423 -0.334 -0.235
mu_bc_H[16] -0.152 0.372 -0.774 -0.194 0.664
mu_bc_R[1] 1.037 0.122 0.800 1.038 1.277
mu_bc_R[2] 1.014 0.075 0.857 1.016 1.155
mu_bc_R[3] 0.220 0.090 0.037 0.224 0.391
mu_bc_R[4] 0.661 0.199 0.257 0.668 1.032
mu_bc_R[5] 1.683 0.377 0.944 1.677 2.427
mu_bc_R[6] -1.444 0.334 -2.094 -1.446 -0.784
mu_bc_R[7] 0.520 0.167 0.180 0.523 0.834
mu_bc_R[8] 0.650 0.149 0.361 0.644 0.948
mu_bc_R[9] 0.384 0.167 0.025 0.390 0.696
mu_bc_R[10] 0.901 0.102 0.709 0.901 1.107
mu_bc_R[11] 1.101 0.074 0.957 1.101 1.247
mu_bc_R[12] 0.864 0.161 0.550 0.862 1.176
mu_bc_R[13] 1.030 0.091 0.841 1.033 1.204
mu_bc_R[14] 0.861 0.122 0.619 0.863 1.096
mu_bc_R[15] 0.780 0.096 0.589 0.781 0.968
mu_bc_R[16] 1.264 0.115 1.032 1.263 1.484
tau_pH[1] 2.448 0.181 2.116 2.442 2.815
tau_pH[2] 1.567 0.148 1.288 1.558 1.875
tau_pH[3] 2.281 0.296 1.749 2.264 2.901
tau_pH[4] 1.113 0.174 0.822 1.097 1.494
tau_pH[5] 0.790 0.128 0.568 0.782 1.070
beta0_pH[1,1] 0.534 0.185 0.133 0.542 0.853
beta0_pH[2,1] 1.424 0.171 1.059 1.437 1.730
beta0_pH[3,1] 1.370 0.243 0.780 1.404 1.750
beta0_pH[4,1] 1.549 0.245 1.028 1.565 1.979
beta0_pH[5,1] -0.785 0.373 -1.687 -0.744 -0.176
beta0_pH[6,1] -0.526 0.439 -1.702 -0.410 0.055
beta0_pH[7,1] 0.753 0.011 0.732 0.753 0.774
beta0_pH[8,1] -0.542 0.301 -1.189 -0.527 -0.025
beta0_pH[9,1] -0.306 0.285 -0.937 -0.275 0.174
beta0_pH[10,1] 0.283 0.189 -0.105 0.292 0.625
beta0_pH[11,1] 0.017 0.319 -0.666 0.030 0.578
beta0_pH[12,1] 0.630 0.193 0.212 0.635 0.985
beta0_pH[13,1] -0.675 0.148 -0.987 -0.669 -0.395
beta0_pH[14,1] 0.007 0.174 -0.356 0.014 0.345
beta0_pH[15,1] 0.063 0.184 -0.331 0.073 0.397
beta0_pH[16,1] -0.874 0.212 -1.305 -0.868 -0.473
beta0_pH[1,2] 2.565 0.086 2.398 2.565 2.737
beta0_pH[2,2] 3.078 0.021 3.037 3.078 3.120
beta0_pH[3,2] 3.626 0.041 3.545 3.625 3.706
beta0_pH[4,2] 2.937 0.019 2.901 2.937 2.974
beta0_pH[5,2] 4.555 1.435 2.755 4.236 8.280
beta0_pH[6,2] 3.333 0.157 3.046 3.331 3.651
beta0_pH[7,2] 2.025 0.049 1.931 2.024 2.124
beta0_pH[8,2] 2.823 0.065 2.698 2.821 2.952
beta0_pH[9,2] 3.966 0.207 3.571 3.962 4.383
beta0_pH[10,2] 4.002 0.130 3.755 3.995 4.261
beta0_pH[11,2] 1.888 0.009 1.870 1.889 1.906
beta0_pH[12,2] 1.801 0.067 1.675 1.799 1.937
beta0_pH[13,2] 2.498 0.026 2.447 2.497 2.548
beta0_pH[14,2] 1.806 0.024 1.760 1.806 1.855
beta0_pH[15,2] 2.550 0.017 2.517 2.550 2.584
beta0_pH[16,2] 2.776 0.025 2.727 2.776 2.824
beta0_pH[1,3] 0.684 0.594 -0.935 0.844 1.369
beta0_pH[2,3] 2.127 0.017 2.094 2.127 2.161
beta0_pH[3,3] 2.689 0.024 2.642 2.689 2.737
beta0_pH[4,3] 2.981 0.034 2.914 2.981 3.048
beta0_pH[5,3] 3.803 2.620 1.007 3.089 11.156
beta0_pH[6,3] 1.260 0.323 0.700 1.233 1.964
beta0_pH[7,3] 0.855 0.049 0.760 0.855 0.951
beta0_pH[8,3] 0.286 0.122 0.030 0.292 0.502
beta0_pH[9,3] -2.916 1.096 -5.302 -2.792 -1.119
beta0_pH[10,3] -0.880 0.891 -3.009 -0.761 0.576
beta0_pH[11,3] -0.003 11.554 -22.520 -0.015 22.520
beta0_pH[12,3] 0.218 11.555 -21.446 0.061 23.470
beta0_pH[13,3] 0.359 11.532 -22.208 0.333 23.960
beta0_pH[14,3] -0.145 11.480 -22.944 -0.177 22.643
beta0_pH[15,3] 0.247 11.636 -22.140 0.079 23.390
beta0_pH[16,3] 0.230 11.433 -22.081 -0.086 23.192
beta0_pH[11,4] 0.991 0.647 -0.507 1.035 2.114
beta0_pH[12,4] 0.788 0.830 -1.172 0.834 2.255
beta0_pH[13,4] 0.730 0.688 -0.941 0.789 1.929
beta0_pH[14,4] 0.325 0.555 -0.977 0.400 1.247
beta0_pH[15,4] 0.634 0.596 -0.789 0.671 1.775
beta0_pH[16,4] 0.956 0.734 -0.647 0.969 2.334
beta0_pH[11,5] 1.669 0.083 1.516 1.665 1.842
beta0_pH[12,5] -1.415 0.459 -2.365 -1.388 -0.547
beta0_pH[13,5] 2.237 0.262 1.810 2.211 2.816
beta0_pH[14,5] 1.879 0.295 1.418 1.843 2.546
beta0_pH[15,5] 1.110 0.065 0.994 1.107 1.243
beta0_pH[16,5] -0.736 0.053 -0.838 -0.737 -0.630
beta1_pH[1,1] 2.628 0.513 1.796 2.565 3.911
beta1_pH[2,1] 1.617 0.327 1.018 1.602 2.323
beta1_pH[3,1] 1.865 0.481 1.133 1.784 3.095
beta1_pH[4,1] 1.941 0.441 1.217 1.889 3.046
beta1_pH[5,1] 2.300 0.487 1.522 2.242 3.495
beta1_pH[6,1] 2.153 0.809 1.106 1.928 4.360
beta1_pH[7,1] 0.000 0.000 0.000 0.000 0.000
beta1_pH[8,1] 2.359 0.812 1.236 2.174 4.369
beta1_pH[9,1] 2.087 0.636 1.233 1.955 3.980
beta1_pH[10,1] 2.144 0.365 1.520 2.127 2.962
beta1_pH[11,1] 3.400 0.516 2.551 3.347 4.571
beta1_pH[12,1] 2.894 0.288 2.334 2.887 3.472
beta1_pH[13,1] 3.484 0.245 3.044 3.475 4.020
beta1_pH[14,1] 2.867 0.317 2.297 2.847 3.509
beta1_pH[15,1] 2.395 0.269 1.893 2.387 2.953
beta1_pH[16,1] 4.041 0.297 3.467 4.042 4.645
beta1_pH[1,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[2,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[3,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[4,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[5,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[6,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[7,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[8,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[9,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[10,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[11,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[12,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[13,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[14,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[15,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[16,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[1,3] 3.459 1.669 1.226 2.993 7.270
beta1_pH[2,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[3,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[4,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[5,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[6,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[7,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[8,3] 4.814 2.334 2.182 4.011 11.178
beta1_pH[9,3] 8.573 2.630 4.402 8.231 14.101
beta1_pH[10,3] 7.170 2.796 3.073 6.618 13.617
beta1_pH[11,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[12,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[13,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[14,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[15,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[16,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[11,4] 2.197 0.697 0.882 2.193 3.677
beta1_pH[12,4] 3.383 0.852 1.837 3.307 5.379
beta1_pH[13,4] 1.725 0.729 0.304 1.706 3.388
beta1_pH[14,4] 2.337 0.632 1.215 2.270 3.784
beta1_pH[15,4] 1.799 0.654 0.475 1.775 3.213
beta1_pH[16,4] 3.908 0.789 2.388 3.883 5.544
beta1_pH[11,5] 0.000 0.000 0.000 0.000 0.000
beta1_pH[12,5] 0.000 0.000 0.000 0.000 0.000
beta1_pH[13,5] 0.000 0.000 0.000 0.000 0.000
beta1_pH[14,5] 0.000 0.000 0.000 0.000 0.000
beta1_pH[15,5] 0.000 0.000 0.000 0.000 0.000
beta1_pH[16,5] 0.000 0.000 0.000 0.000 0.000
beta2_pH[1,1] 1.150 0.808 0.303 0.862 3.332
beta2_pH[2,1] 1.557 0.938 0.289 1.335 3.407
beta2_pH[3,1] 1.269 0.832 0.196 1.098 2.903
beta2_pH[4,1] 1.089 0.762 0.217 0.828 2.877
beta2_pH[5,1] 2.198 1.792 0.286 1.754 7.344
beta2_pH[6,1] 0.995 1.440 0.114 0.427 5.483
beta2_pH[7,1] 0.000 0.000 0.000 0.000 0.000
beta2_pH[8,1] 0.727 0.953 0.139 0.367 3.802
beta2_pH[9,1] 0.973 1.367 0.149 0.509 5.189
beta2_pH[10,1] 1.273 1.276 0.292 0.851 5.202
beta2_pH[11,1] 0.409 0.230 0.184 0.350 0.958
beta2_pH[12,1] 1.547 0.794 0.625 1.287 3.628
beta2_pH[13,1] 0.978 0.385 0.522 0.892 1.972
beta2_pH[14,1] 1.105 0.645 0.433 0.912 3.214
beta2_pH[15,1] 1.052 0.606 0.398 0.879 2.847
beta2_pH[16,1] 1.046 0.346 0.622 0.970 1.874
beta2_pH[1,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[2,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[3,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[4,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[5,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[6,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[7,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[8,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[9,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[10,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[11,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[12,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[13,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[14,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[15,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[16,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[1,3] 2.251 4.014 0.114 0.493 14.928
beta2_pH[2,3] 0.000 0.000 0.000 0.000 0.000
beta2_pH[3,3] 0.000 0.000 0.000 0.000 0.000
beta2_pH[4,3] 0.000 0.000 0.000 0.000 0.000
beta2_pH[5,3] 0.000 0.000 0.000 0.000 0.000
beta2_pH[6,3] 0.000 0.000 0.000 0.000 0.000
beta2_pH[7,3] 0.000 0.000 0.000 0.000 0.000
beta2_pH[8,3] 2.797 2.305 0.551 1.840 8.802
beta2_pH[9,3] 0.178 0.100 0.080 0.155 0.388
beta2_pH[10,3] 0.329 0.468 0.099 0.232 0.999
beta2_pH[11,3] -8.277 5.870 -22.204 -7.169 -0.442
beta2_pH[12,3] -8.223 5.833 -22.301 -7.179 -0.421
beta2_pH[13,3] -8.381 5.912 -22.286 -7.204 -0.370
beta2_pH[14,3] -8.385 5.897 -22.054 -7.222 -0.391
beta2_pH[15,3] -8.314 5.883 -22.382 -7.212 -0.405
beta2_pH[16,3] -8.329 5.902 -22.296 -7.232 -0.404
beta2_pH[11,4] -4.357 3.591 -14.546 -3.381 -0.638
beta2_pH[12,4] -4.917 3.617 -14.520 -3.916 -1.023
beta2_pH[13,4] -4.507 3.630 -14.779 -3.455 -0.779
beta2_pH[14,4] -4.556 3.488 -14.529 -3.533 -0.971
beta2_pH[15,4] -4.187 3.357 -14.037 -3.299 -0.754
beta2_pH[16,4] -4.655 3.353 -14.399 -3.633 -1.281
beta2_pH[11,5] 0.000 0.000 0.000 0.000 0.000
beta2_pH[12,5] 0.000 0.000 0.000 0.000 0.000
beta2_pH[13,5] 0.000 0.000 0.000 0.000 0.000
beta2_pH[14,5] 0.000 0.000 0.000 0.000 0.000
beta2_pH[15,5] 0.000 0.000 0.000 0.000 0.000
beta2_pH[16,5] 0.000 0.000 0.000 0.000 0.000
beta3_pH[1,1] 34.980 1.137 32.878 34.973 37.469
beta3_pH[2,1] 32.890 1.504 30.368 32.697 36.382
beta3_pH[3,1] 33.052 1.705 29.781 32.955 37.020
beta3_pH[4,1] 32.633 1.894 29.589 32.451 37.024
beta3_pH[5,1] 27.604 1.079 26.043 27.499 30.327
beta3_pH[6,1] 33.247 3.171 26.657 33.108 40.810
beta3_pH[7,1] 25.000 0.000 25.000 25.000 25.000
beta3_pH[8,1] 35.882 3.546 29.563 36.128 42.022
beta3_pH[9,1] 32.608 2.915 27.845 32.256 40.086
beta3_pH[10,1] 32.627 1.351 30.218 32.525 35.600
beta3_pH[11,1] 31.158 1.145 29.246 31.112 33.599
beta3_pH[12,1] 30.763 0.482 29.794 30.766 31.663
beta3_pH[13,1] 31.400 0.530 30.407 31.389 32.467
beta3_pH[14,1] 32.545 0.757 31.266 32.487 34.208
beta3_pH[15,1] 30.968 0.703 29.571 30.963 32.332
beta3_pH[16,1] 30.128 0.347 29.470 30.124 30.800
beta3_pH[1,2] 34.391 5.456 25.471 34.234 43.504
beta3_pH[2,2] 34.349 5.471 25.496 34.174 43.473
beta3_pH[3,2] 30.405 7.156 18.663 30.333 42.314
beta3_pH[4,2] 30.692 7.257 18.638 30.707 42.390
beta3_pH[5,2] 30.407 7.097 18.730 30.344 42.322
beta3_pH[6,2] 30.528 7.218 18.606 30.562 42.464
beta3_pH[7,2] 25.000 0.000 25.000 25.000 25.000
beta3_pH[8,2] 25.000 0.000 25.000 25.000 25.000
beta3_pH[9,2] 25.000 0.000 25.000 25.000 25.000
beta3_pH[10,2] 25.000 0.000 25.000 25.000 25.000
beta3_pH[11,2] 24.000 0.000 24.000 24.000 24.000
beta3_pH[12,2] 24.000 0.000 24.000 24.000 24.000
beta3_pH[13,2] 24.000 0.000 24.000 24.000 24.000
beta3_pH[14,2] 24.000 0.000 24.000 24.000 24.000
beta3_pH[15,2] 24.000 0.000 24.000 24.000 24.000
beta3_pH[16,2] 24.000 0.000 24.000 24.000 24.000
beta3_pH[1,3] 39.781 2.596 33.282 40.157 43.666
beta3_pH[2,3] 34.447 5.511 25.447 34.405 43.444
beta3_pH[3,3] 30.675 7.245 18.591 30.864 42.436
beta3_pH[4,3] 30.379 7.280 18.506 30.381 42.371
beta3_pH[5,3] 25.000 0.000 25.000 25.000 25.000
beta3_pH[6,3] 25.000 0.000 25.000 25.000 25.000
beta3_pH[7,3] 25.000 0.000 25.000 25.000 25.000
beta3_pH[8,3] 42.112 0.834 41.010 41.817 43.887
beta3_pH[9,3] 36.000 1.749 34.528 35.347 41.299
beta3_pH[10,3] 37.774 2.544 34.619 37.035 43.330
beta3_pH[11,3] 38.412 2.572 34.247 38.392 42.742
beta3_pH[12,3] 38.560 2.603 34.262 38.578 42.775
beta3_pH[13,3] 38.437 2.618 34.195 38.415 42.813
beta3_pH[14,3] 38.530 2.611 34.224 38.531 42.787
beta3_pH[15,3] 38.479 2.601 34.242 38.471 42.741
beta3_pH[16,3] 38.533 2.630 34.184 38.650 42.790
beta3_pH[11,4] 40.724 1.139 38.738 40.597 43.219
beta3_pH[12,4] 41.934 0.763 40.281 42.168 42.942
beta3_pH[13,4] 41.529 1.133 39.232 41.490 43.610
beta3_pH[14,4] 40.339 0.744 38.834 40.352 41.713
beta3_pH[15,4] 41.182 0.984 39.206 41.103 43.302
beta3_pH[16,4] 41.758 0.410 40.941 41.799 42.555
beta3_pH[11,5] 25.000 0.000 25.000 25.000 25.000
beta3_pH[12,5] 25.000 0.000 25.000 25.000 25.000
beta3_pH[13,5] 25.000 0.000 25.000 25.000 25.000
beta3_pH[14,5] 25.000 0.000 25.000 25.000 25.000
beta3_pH[15,5] 25.000 0.000 25.000 25.000 25.000
beta3_pH[16,5] 25.000 0.000 25.000 25.000 25.000
beta4_pH[1,1] 0.262 0.021 0.221 0.263 0.303
beta4_pH[2,1] 0.211 0.005 0.200 0.211 0.220
beta4_pH[3,1] 0.173 0.007 0.159 0.173 0.186
beta4_pH[4,1] 0.294 0.011 0.272 0.293 0.316
beta4_pH[5,1] 0.427 0.251 0.092 0.394 0.954
beta4_pH[6,1] 0.776 0.389 0.334 0.733 1.402
beta4_pH[7,1] 0.855 0.345 0.354 0.798 1.691
beta4_pH[8,1] 0.159 0.070 0.045 0.151 0.314
beta4_pH[9,1] 0.214 0.088 0.062 0.207 0.400
beta4_pH[10,1] 0.001 0.001 0.000 0.001 0.003
beta4_pH[11,1] 0.446 0.012 0.424 0.446 0.471
beta4_pH[12,1] 0.642 0.022 0.599 0.642 0.687
beta4_pH[13,1] 0.630 0.014 0.603 0.630 0.657
beta4_pH[14,1] 0.667 0.031 0.608 0.667 0.727
beta4_pH[15,1] 0.796 0.021 0.756 0.796 0.837
beta4_pH[16,1] 1.216 0.042 1.137 1.215 1.303
beta4_pH[1,2] 0.295 0.095 0.127 0.289 0.501
beta4_pH[2,2] 0.515 0.022 0.471 0.514 0.560
beta4_pH[3,2] 0.446 0.012 0.422 0.446 0.469
beta4_pH[4,2] 0.653 0.038 0.581 0.654 0.726
beta4_pH[5,2] 0.700 0.685 0.109 0.525 2.272
beta4_pH[6,2] 0.891 0.758 0.210 0.733 2.651
beta4_pH[7,2] 0.894 0.456 0.343 0.815 1.945
beta4_pH[8,2] 0.854 0.202 0.489 0.844 1.274
beta4_pH[9,2] 1.165 1.194 0.316 0.869 3.800
beta4_pH[10,2] 0.458 0.052 0.358 0.457 0.563
beta4_pH[11,2] 1.268 0.050 1.171 1.268 1.368
beta4_pH[12,2] 0.603 0.112 0.395 0.600 0.834
beta4_pH[13,2] 0.948 0.045 0.862 0.948 1.037
beta4_pH[14,2] 0.931 0.062 0.813 0.929 1.057
beta4_pH[15,2] 1.010 0.025 0.961 1.010 1.060
beta4_pH[16,2] 1.392 0.075 1.248 1.390 1.545
beta4_pH[1,3] 0.239 0.075 0.104 0.235 0.398
beta4_pH[2,3] 0.050 0.017 0.020 0.049 0.085
beta4_pH[3,3] 0.052 0.010 0.033 0.052 0.071
beta4_pH[4,3] 0.167 0.031 0.108 0.166 0.226
beta4_pH[5,3] 0.766 0.883 0.115 0.561 2.786
beta4_pH[6,3] 0.898 0.767 0.214 0.727 2.573
beta4_pH[7,3] 0.904 0.450 0.364 0.829 1.874
beta4_pH[8,3] 0.393 0.184 0.103 0.370 0.824
beta4_pH[9,3] 0.906 0.740 0.155 0.713 2.815
beta4_pH[10,3] 0.055 0.037 0.003 0.048 0.140
beta4_pH[11,3] 0.616 0.258 0.264 0.570 1.245
beta4_pH[12,3] 0.703 0.225 0.369 0.670 1.258
beta4_pH[13,3] 0.713 0.186 0.414 0.690 1.121
beta4_pH[14,3] 0.859 0.223 0.510 0.827 1.390
beta4_pH[15,3] 0.824 0.254 0.433 0.784 1.422
beta4_pH[16,3] 0.833 0.383 0.333 0.762 1.765
beta4_pH[1,4] 0.647 0.403 0.198 0.565 1.588
beta4_pH[2,4] 0.821 0.879 0.111 0.584 2.896
beta4_pH[3,4] 0.884 1.015 0.116 0.616 3.411
beta4_pH[4,4] 0.737 0.467 0.205 0.627 1.861
beta4_pH[5,4] 0.761 0.818 0.119 0.553 2.674
beta4_pH[6,4] 0.913 0.781 0.216 0.735 2.629
beta4_pH[7,4] 0.894 0.421 0.345 0.805 1.914
beta4_pH[8,4] 0.750 0.720 0.152 0.584 2.397
beta4_pH[9,4] 0.775 1.004 0.104 0.560 2.686
beta4_pH[10,4] 1.458 3.146 0.006 0.370 11.271
beta4_pH[11,4] 0.395 0.016 0.365 0.395 0.427
beta4_pH[12,4] 0.550 0.054 0.446 0.549 0.655
beta4_pH[13,4] 0.835 0.040 0.761 0.834 0.917
beta4_pH[14,4] 0.864 0.072 0.728 0.862 1.009
beta4_pH[15,4] 0.626 0.017 0.593 0.626 0.660
beta4_pH[16,4] 0.647 0.039 0.575 0.646 0.726
beta4_pH[1,5] 0.636 0.382 0.204 0.554 1.630
beta4_pH[2,5] 0.822 0.857 0.119 0.608 2.811
beta4_pH[3,5] 0.881 0.996 0.114 0.605 3.341
beta4_pH[4,5] 0.733 0.484 0.211 0.628 1.879
beta4_pH[5,5] 0.783 1.046 0.112 0.560 2.806
beta4_pH[6,5] 0.918 0.825 0.212 0.736 2.672
beta4_pH[7,5] 0.879 0.388 0.343 0.805 1.882
beta4_pH[8,5] 0.739 0.748 0.151 0.584 2.117
beta4_pH[9,5] 0.763 0.728 0.124 0.571 2.690
beta4_pH[10,5] 1.542 3.202 0.009 0.399 12.089
beta4_pH[11,5] 0.362 0.075 0.214 0.362 0.513
beta4_pH[12,5] 0.703 0.210 0.381 0.672 1.201
beta4_pH[13,5] 0.970 0.111 0.768 0.964 1.201
beta4_pH[14,5] 0.845 0.156 0.574 0.833 1.178
beta4_pH[15,5] 1.305 0.116 1.080 1.302 1.540
beta4_pH[16,5] 0.927 0.380 0.372 0.861 1.860
beta0_pelagic[1] 2.268 0.036 2.200 2.266 2.340
beta0_pelagic[2] 1.537 0.007 1.524 1.537 1.551
beta0_pelagic[3] 0.202 0.249 -0.449 0.257 0.550
beta0_pelagic[4] 0.224 0.232 -0.324 0.236 0.669
beta0_pelagic[5] 1.444 0.153 1.157 1.442 1.747
beta0_pelagic[6] 1.355 0.060 1.244 1.350 1.481
beta0_pelagic[7] 1.471 0.038 1.401 1.469 1.549
beta0_pelagic[8] 1.671 0.044 1.588 1.669 1.758
beta0_pelagic[9] 2.776 0.063 2.660 2.774 2.903
beta0_pelagic[10] 2.603 0.032 2.544 2.602 2.667
beta0_pelagic[11] -0.168 0.205 -0.638 -0.146 0.169
beta0_pelagic[12] 1.801 0.037 1.734 1.799 1.875
beta0_pelagic[13] -1.616 0.780 -3.306 -1.498 -0.381
beta0_pelagic[14] -0.275 0.168 -0.665 -0.258 0.012
beta0_pelagic[15] -0.209 0.033 -0.299 -0.195 -0.180
beta0_pelagic[16] -0.482 0.377 -1.333 -0.390 0.067
beta1_pelagic[1] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[2] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[3] 0.762 0.289 0.339 0.717 1.645
beta1_pelagic[4] 0.825 0.259 0.330 0.816 1.430
beta1_pelagic[5] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[6] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[7] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[8] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[9] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[10] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[11] 1.189 0.284 0.719 1.155 1.856
beta1_pelagic[12] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[13] 3.981 1.479 1.975 3.588 7.127
beta1_pelagic[14] 1.779 0.447 1.124 1.683 2.946
beta1_pelagic[15] 0.091 0.155 0.000 0.014 0.579
beta1_pelagic[16] 1.620 0.572 0.774 1.507 2.822
beta2_pelagic[1] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[2] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[3] 6.342 6.018 0.174 4.607 21.640
beta2_pelagic[4] 7.488 6.131 0.332 6.116 21.856
beta2_pelagic[5] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[6] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[7] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[8] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[9] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[10] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[11] 1.704 1.494 0.256 1.169 5.754
beta2_pelagic[12] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[13] 0.133 0.061 0.056 0.123 0.281
beta2_pelagic[14] 0.875 0.869 0.201 0.603 3.581
beta2_pelagic[15] 2.629 1.883 0.175 2.284 7.134
beta2_pelagic[16] 0.716 0.968 0.120 0.353 3.956
beta3_pelagic[1] 28.750 3.912 22.326 28.648 35.213
beta3_pelagic[2] 28.863 3.928 22.313 28.981 35.195
beta3_pelagic[3] 28.163 2.705 23.327 28.702 32.811
beta3_pelagic[4] 26.151 1.483 23.165 26.276 29.644
beta3_pelagic[5] 29.692 6.630 18.759 29.715 40.468
beta3_pelagic[6] 29.670 6.653 18.699 29.587 40.495
beta3_pelagic[7] 29.383 6.576 18.555 29.195 40.403
beta3_pelagic[8] 29.624 6.758 18.602 29.614 40.562
beta3_pelagic[9] 29.553 6.569 18.634 29.452 40.349
beta3_pelagic[10] 29.361 6.637 18.509 29.122 40.428
beta3_pelagic[11] 30.069 1.296 27.704 29.992 32.987
beta3_pelagic[12] 30.000 0.000 30.000 30.000 30.000
beta3_pelagic[13] 30.989 4.502 23.631 30.652 40.447
beta3_pelagic[14] 34.245 1.708 31.471 34.011 38.366
beta3_pelagic[15] 34.271 5.770 23.608 35.217 41.760
beta3_pelagic[16] 31.275 2.655 26.329 31.175 37.213
beta0_yellow[1] -0.484 0.080 -0.637 -0.485 -0.329
beta0_yellow[2] 0.540 0.085 0.352 0.548 0.686
beta0_yellow[3] -0.208 0.131 -0.465 -0.211 0.013
beta0_yellow[4] 0.904 0.246 0.212 0.975 1.193
beta0_yellow[5] 0.246 0.327 -0.458 0.270 0.823
beta0_yellow[6] 0.822 0.105 0.620 0.822 1.026
beta0_yellow[7] 1.013 0.102 0.812 1.011 1.220
beta0_yellow[8] 1.017 0.095 0.834 1.016 1.207
beta0_yellow[9] 0.442 0.098 0.242 0.445 0.633
beta0_yellow[10] 0.540 0.070 0.402 0.541 0.675
beta0_yellow[11] -0.192 0.121 -0.457 -0.184 0.014
beta0_yellow[12] -1.085 0.094 -1.276 -1.085 -0.903
beta0_yellow[13] -1.338 0.169 -1.721 -1.321 -1.061
beta0_yellow[14] -0.249 0.154 -0.580 -0.231 -0.007
beta0_yellow[15] -1.243 0.132 -1.561 -1.229 -1.031
beta0_yellow[16] -0.256 0.091 -0.436 -0.255 -0.077
beta1_yellow[1] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[2] 1.017 0.236 0.575 1.012 1.479
beta1_yellow[3] 0.466 0.248 0.019 0.476 0.922
beta1_yellow[4] 1.310 0.562 0.642 1.181 3.024
beta1_yellow[5] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[6] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[7] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[8] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[9] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[10] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[11] 0.251 0.335 0.001 0.066 1.105
beta1_yellow[12] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[13] 1.915 1.680 0.056 1.135 5.938
beta1_yellow[14] 0.740 1.019 0.012 0.433 3.991
beta1_yellow[15] 0.159 0.195 0.011 0.075 0.701
beta1_yellow[16] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[1] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[2] -7.563 5.736 -21.575 -6.376 -0.428
beta2_yellow[3] -7.829 5.821 -22.218 -6.735 -0.372
beta2_yellow[4] -5.923 6.055 -21.078 -4.081 -0.126
beta2_yellow[5] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[6] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[7] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[8] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[9] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[10] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[11] -7.798 5.708 -21.055 -6.618 -0.392
beta2_yellow[12] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[13] -5.294 5.303 -19.365 -3.271 -0.313
beta2_yellow[14] -7.227 5.729 -20.568 -6.035 -0.365
beta2_yellow[15] -7.577 5.654 -21.001 -6.397 -0.334
beta2_yellow[16] 0.000 0.000 0.000 0.000 0.000
beta3_yellow[1] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[2] 28.888 1.150 27.325 28.674 31.853
beta3_yellow[3] 32.307 3.286 22.949 32.494 38.476
beta3_yellow[4] 29.852 3.204 24.400 28.617 35.612
beta3_yellow[5] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[6] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[7] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[8] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[9] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[10] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[11] 35.332 3.774 30.381 34.328 43.342
beta3_yellow[12] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[13] 32.347 2.696 29.087 32.049 37.969
beta3_yellow[14] 33.808 3.523 29.178 33.385 42.751
beta3_yellow[15] 36.521 3.818 30.273 36.303 43.462
beta3_yellow[16] 35.000 0.000 35.000 35.000 35.000
beta0_black[1] 0.009 0.074 -0.139 0.009 0.155
beta0_black[2] 1.856 0.032 1.793 1.856 1.917
beta0_black[3] 1.240 0.038 1.165 1.239 1.313
beta0_black[4] 2.425 0.047 2.335 2.423 2.518
beta0_black[5] 5.151 2.393 1.801 4.688 11.214
beta0_black[6] 5.221 2.120 2.634 4.712 10.690
beta0_black[7] 4.148 2.197 1.635 3.548 10.291
beta0_black[8] 0.886 0.213 0.495 0.877 1.335
beta0_black[9] 2.802 0.271 2.275 2.801 3.345
beta0_black[10] 1.288 0.050 1.191 1.287 1.388
beta0_black[11] 2.977 0.052 2.876 2.977 3.081
beta0_black[12] 4.910 0.190 4.549 4.914 5.282
beta0_black[13] 0.449 0.339 0.016 0.371 1.242
beta0_black[14] 2.851 0.108 2.641 2.850 3.067
beta0_black[15] 0.723 0.052 0.624 0.723 0.824
beta0_black[16] 3.774 0.074 3.634 3.773 3.922
beta2_black[1] 0.000 0.000 0.000 0.000 0.000
beta2_black[2] 0.000 0.000 0.000 0.000 0.000
beta2_black[3] 0.000 0.000 0.000 0.000 0.000
beta2_black[4] 0.000 0.000 0.000 0.000 0.000
beta2_black[5] 0.000 0.000 0.000 0.000 0.000
beta2_black[6] 0.000 0.000 0.000 0.000 0.000
beta2_black[7] 0.000 0.000 0.000 0.000 0.000
beta2_black[8] 0.000 0.000 0.000 0.000 0.000
beta2_black[9] 0.000 0.000 0.000 0.000 0.000
beta2_black[10] 0.000 0.000 0.000 0.000 0.000
beta2_black[11] 0.000 0.000 0.000 0.000 0.000
beta2_black[12] 0.000 0.000 0.000 0.000 0.000
beta2_black[13] -2.589 1.801 -7.155 -2.179 -0.308
beta2_black[14] 0.000 0.000 0.000 0.000 0.000
beta2_black[15] 0.000 0.000 0.000 0.000 0.000
beta2_black[16] 0.000 0.000 0.000 0.000 0.000
beta3_black[1] 25.000 0.000 25.000 25.000 25.000
beta3_black[2] 25.000 0.000 25.000 25.000 25.000
beta3_black[3] 25.000 0.000 25.000 25.000 25.000
beta3_black[4] 25.000 0.000 25.000 25.000 25.000
beta3_black[5] 25.000 0.000 25.000 25.000 25.000
beta3_black[6] 25.000 0.000 25.000 25.000 25.000
beta3_black[7] 25.000 0.000 25.000 25.000 25.000
beta3_black[8] 25.000 0.000 25.000 25.000 25.000
beta3_black[9] 25.000 0.000 25.000 25.000 25.000
beta3_black[10] 25.000 0.000 25.000 25.000 25.000
beta3_black[11] 25.000 0.000 25.000 25.000 25.000
beta3_black[12] 25.000 0.000 25.000 25.000 25.000
beta3_black[13] 37.299 3.413 27.189 38.094 40.841
beta3_black[14] 25.000 0.000 25.000 25.000 25.000
beta3_black[15] 25.000 0.000 25.000 25.000 25.000
beta3_black[16] 25.000 0.000 25.000 25.000 25.000
beta4_black[1] -0.324 0.091 -0.499 -0.324 -0.141
beta4_black[2] 0.257 0.052 0.154 0.258 0.358
beta4_black[3] -0.721 0.064 -0.848 -0.720 -0.596
beta4_black[4] 0.629 0.130 0.382 0.627 0.888
beta4_black[5] 0.745 1.452 -1.139 0.495 4.115
beta4_black[6] 0.824 1.422 -0.904 0.536 4.228
beta4_black[7] 0.775 1.590 -1.052 0.509 4.060
beta4_black[8] 0.073 0.285 -0.520 0.082 0.615
beta4_black[9] 1.024 0.784 -0.052 0.873 2.941
beta4_black[10] 0.072 0.067 -0.062 0.073 0.204
beta4_black[11] -0.610 0.065 -0.736 -0.609 -0.484
beta4_black[12] 1.248 0.286 0.710 1.242 1.807
beta4_black[13] -1.003 0.076 -1.159 -1.002 -0.858
beta4_black[14] 0.229 0.185 -0.120 0.226 0.604
beta4_black[15] 0.001 0.046 -0.088 0.002 0.090
beta4_black[16] -0.420 0.115 -0.649 -0.419 -0.196
beta0_dsr[11] -0.337 0.009 -0.355 -0.337 -0.318
beta0_dsr[12] 1.422 0.043 1.341 1.421 1.507
beta0_dsr[13] 0.131 0.022 0.087 0.130 0.177
beta0_dsr[14] -0.135 0.017 -0.170 -0.135 -0.101
beta0_dsr[15] 0.185 0.014 0.157 0.184 0.212
beta0_dsr[16] 0.208 0.015 0.178 0.208 0.238
beta1_dsr[11] 0.000 0.000 0.000 0.000 0.000
beta1_dsr[12] 0.000 0.000 0.000 0.000 0.000
beta1_dsr[13] 0.000 0.000 0.000 0.000 0.000
beta1_dsr[14] 0.000 0.000 0.000 0.000 0.000
beta1_dsr[15] 0.000 0.000 0.000 0.000 0.000
beta1_dsr[16] 0.000 0.000 0.000 0.000 0.000
beta4_dsr[11] 0.548 0.036 0.477 0.548 0.615
beta4_dsr[12] 0.811 0.240 0.353 0.805 1.300
beta4_dsr[13] -0.046 0.042 -0.130 -0.046 0.035
beta4_dsr[14] -0.206 0.072 -0.347 -0.207 -0.065
beta4_dsr[15] 0.236 0.022 0.193 0.236 0.279
beta4_dsr[16] -0.098 0.045 -0.186 -0.098 -0.009
beta0_slope[11] -2.200 0.021 -2.242 -2.200 -2.160
beta0_slope[12] -5.431 0.231 -5.879 -5.430 -4.985
beta0_slope[13] -1.455 0.034 -1.522 -1.454 -1.390
beta0_slope[14] -2.692 0.051 -2.793 -2.693 -2.594
beta0_slope[15] -1.567 0.021 -1.609 -1.566 -1.524
beta0_slope[16] -3.014 0.046 -3.103 -3.013 -2.925
beta1_slope[11] 0.000 0.000 0.000 0.000 0.000
beta1_slope[12] 0.000 0.000 0.000 0.000 0.000
beta1_slope[13] 0.000 0.000 0.000 0.000 0.000
beta1_slope[14] 0.000 0.000 0.000 0.000 0.000
beta1_slope[15] 0.000 0.000 0.000 0.000 0.000
beta1_slope[16] 0.000 0.000 0.000 0.000 0.000
beta4_slope[11] -0.343 0.075 -0.498 -0.341 -0.205
beta4_slope[12] -4.673 0.798 -6.352 -4.647 -3.235
beta4_slope[13] 0.280 0.053 0.179 0.278 0.386
beta4_slope[14] -0.178 0.152 -0.492 -0.172 0.109
beta4_slope[15] -0.386 0.033 -0.451 -0.386 -0.319
beta4_slope[16] 0.219 0.106 0.006 0.223 0.429
beta5_slope[11] 3.078 0.107 2.882 3.076 3.298
beta5_slope[12] 2.017 0.383 1.344 2.000 2.825
beta5_slope[13] 0.935 0.079 0.779 0.934 1.093
beta5_slope[14] 3.792 0.359 3.114 3.781 4.544
beta5_slope[15] 1.387 0.074 1.243 1.388 1.536
beta5_slope[16] 2.694 0.233 2.257 2.693 3.163
tau_comp[1] 3.768 0.287 3.235 3.756 4.350
tau_comp[2] 2.846 0.249 2.397 2.834 3.356
tau_comp[3] 0.868 0.073 0.735 0.863 1.020
tau_comp[4] 1.270 0.195 0.939 1.254 1.692
tau_comp[5] 1.202 0.179 0.878 1.194 1.581
sigma_H[1] 0.205 0.049 0.117 0.202 0.308
sigma_H[2] 0.174 0.028 0.127 0.171 0.237
sigma_H[3] 0.188 0.040 0.115 0.186 0.276
sigma_H[4] 0.465 0.081 0.332 0.457 0.650
sigma_H[5] 0.963 0.201 0.601 0.953 1.386
sigma_H[6] 0.528 0.196 0.156 0.517 0.931
sigma_H[7] 0.309 0.061 0.213 0.300 0.453
sigma_H[8] 0.424 0.096 0.265 0.415 0.625
sigma_H[9] 0.517 0.119 0.331 0.499 0.790
sigma_H[10] 0.209 0.038 0.143 0.206 0.292
sigma_H[11] 0.266 0.043 0.192 0.263 0.364
sigma_H[12] 0.381 0.145 0.199 0.336 0.713
sigma_H[13] 0.217 0.036 0.157 0.213 0.301
sigma_H[14] 0.497 0.090 0.343 0.489 0.692
sigma_H[15] 0.235 0.038 0.172 0.231 0.321
sigma_H[16] 0.212 0.041 0.146 0.208 0.305
lambda_H[1] 2.206 3.215 0.076 1.036 11.558
lambda_H[2] 5.180 5.040 0.495 3.714 18.318
lambda_H[3] 5.405 8.096 0.294 2.782 25.723
lambda_H[4] 0.006 0.004 0.001 0.005 0.017
lambda_H[5] 2.750 7.444 0.028 0.698 20.573
lambda_H[6] 6.240 13.571 0.009 0.254 45.510
lambda_H[7] 0.012 0.008 0.002 0.010 0.033
lambda_H[8] 7.538 9.932 0.045 4.194 35.229
lambda_H[9] 0.015 0.010 0.003 0.013 0.041
lambda_H[10] 0.220 0.306 0.025 0.151 0.797
lambda_H[11] 0.227 0.346 0.011 0.106 1.117
lambda_H[12] 3.715 5.088 0.179 2.042 18.067
lambda_H[13] 3.148 2.971 0.260 2.283 10.747
lambda_H[14] 3.238 3.922 0.200 2.121 14.107
lambda_H[15] 0.021 0.099 0.003 0.014 0.065
lambda_H[16] 0.501 0.789 0.032 0.258 2.424
mu_lambda_H[1] 4.071 1.919 1.105 3.833 8.259
mu_lambda_H[2] 3.682 1.977 0.528 3.493 7.817
mu_lambda_H[3] 3.339 1.860 0.642 3.037 7.665
sigma_lambda_H[1] 8.148 4.382 1.782 7.400 18.018
sigma_lambda_H[2] 8.150 4.783 0.874 7.609 18.347
sigma_lambda_H[3] 6.166 4.089 0.894 5.249 16.388
beta_H[1,1] 6.483 1.354 3.149 6.710 8.435
beta_H[2,1] 9.747 0.626 8.334 9.785 10.889
beta_H[3,1] 7.819 0.801 5.882 7.926 9.084
beta_H[4,1] 10.495 7.725 -5.149 10.751 25.708
beta_H[5,1] -0.278 2.549 -5.508 -0.122 4.184
beta_H[6,1] 3.008 4.370 -7.776 4.449 8.073
beta_H[7,1] 1.265 5.642 -10.402 1.568 11.652
beta_H[8,1] 1.460 4.229 -2.651 1.228 4.449
beta_H[9,1] 13.546 5.611 2.590 13.458 24.714
beta_H[10,1] 6.851 1.849 2.865 6.877 10.392
beta_H[11,1] 4.824 3.518 -2.862 5.331 9.875
beta_H[12,1] 2.614 1.113 0.677 2.495 5.291
beta_H[13,1] 9.087 0.867 7.305 9.130 10.546
beta_H[14,1] 2.131 1.055 0.168 2.139 4.108
beta_H[15,1] -6.960 3.506 -13.422 -7.066 0.718
beta_H[16,1] 3.930 3.029 -0.785 3.489 11.091
beta_H[1,2] 7.784 0.271 7.252 7.792 8.296
beta_H[2,2] 9.881 0.163 9.558 9.883 10.191
beta_H[3,2] 8.798 0.203 8.401 8.794 9.202
beta_H[4,2] 3.027 1.456 0.277 2.980 6.109
beta_H[5,2] 1.733 1.026 -0.293 1.726 3.717
beta_H[6,2] 6.025 1.098 3.536 6.190 7.768
beta_H[7,2] 2.378 1.048 0.446 2.346 4.477
beta_H[8,2] 2.968 1.202 0.998 3.140 4.275
beta_H[9,2] 3.199 1.080 1.128 3.174 5.364
beta_H[10,2] 8.033 0.367 7.285 8.038 8.750
beta_H[11,2] 9.832 0.636 8.863 9.757 11.207
beta_H[12,2] 3.937 0.372 3.231 3.920 4.697
beta_H[13,2] 9.157 0.241 8.712 9.152 9.651
beta_H[14,2] 4.019 0.340 3.326 4.015 4.686
beta_H[15,2] 11.528 0.635 10.185 11.532 12.708
beta_H[16,2] 4.478 0.826 2.983 4.465 6.166
beta_H[1,3] 8.443 0.259 7.987 8.434 8.969
beta_H[2,3] 9.861 0.145 9.587 9.856 10.157
beta_H[3,3] 9.460 0.163 9.142 9.457 9.793
beta_H[4,3] -2.187 0.862 -3.930 -2.180 -0.521
beta_H[5,3] 3.894 0.614 2.599 3.894 5.079
beta_H[6,3] 8.404 1.189 6.721 8.100 10.855
beta_H[7,3] -2.553 0.622 -3.821 -2.557 -1.360
beta_H[8,3] 5.353 0.580 4.715 5.260 6.702
beta_H[9,3] -2.725 0.699 -4.103 -2.717 -1.393
beta_H[10,3] 8.550 0.279 8.010 8.552 9.099
beta_H[11,3] 8.499 0.285 7.919 8.521 9.000
beta_H[12,3] 5.268 0.323 4.509 5.306 5.793
beta_H[13,3] 8.861 0.179 8.493 8.865 9.206
beta_H[14,3] 5.778 0.275 5.173 5.802 6.257
beta_H[15,3] 10.281 0.299 9.704 10.283 10.850
beta_H[16,3] 5.922 0.585 4.748 5.937 7.012
beta_H[1,4] 8.098 0.188 7.702 8.111 8.434
beta_H[2,4] 9.849 0.130 9.576 9.851 10.097
beta_H[3,4] 9.949 0.150 9.630 9.961 10.219
beta_H[4,4] 11.814 0.462 10.939 11.814 12.731
beta_H[5,4] 5.569 0.730 4.352 5.509 7.162
beta_H[6,4] 7.189 0.927 5.237 7.405 8.479
beta_H[7,4] 8.542 0.345 7.862 8.544 9.200
beta_H[8,4] 6.726 0.284 6.066 6.753 7.177
beta_H[9,4] 7.437 0.449 6.564 7.436 8.320
beta_H[10,4] 7.523 0.219 7.115 7.518 7.965
beta_H[11,4] 9.500 0.192 9.131 9.498 9.887
beta_H[12,4] 7.261 0.202 6.881 7.254 7.691
beta_H[13,4] 9.058 0.138 8.778 9.066 9.317
beta_H[14,4] 7.828 0.214 7.415 7.825 8.267
beta_H[15,4] 9.474 0.222 9.037 9.473 9.895
beta_H[16,4] 9.431 0.236 8.991 9.427 9.913
beta_H[1,5] 8.900 0.152 8.582 8.903 9.200
beta_H[2,5] 10.676 0.100 10.489 10.672 10.888
beta_H[3,5] 10.761 0.161 10.458 10.756 11.091
beta_H[4,5] 8.349 0.486 7.437 8.335 9.315
beta_H[5,5] 5.203 0.610 3.723 5.281 6.166
beta_H[6,5] 8.991 0.633 8.024 8.871 10.399
beta_H[7,5] 6.571 0.343 5.898 6.574 7.236
beta_H[8,5] 8.261 0.234 7.901 8.235 8.798
beta_H[9,5] 8.174 0.458 7.258 8.173 9.073
beta_H[10,5] 10.045 0.224 9.591 10.050 10.485
beta_H[11,5] 11.594 0.219 11.164 11.594 12.027
beta_H[12,5] 8.607 0.190 8.233 8.606 8.978
beta_H[13,5] 10.134 0.132 9.891 10.130 10.402
beta_H[14,5] 9.235 0.211 8.841 9.233 9.684
beta_H[15,5] 11.307 0.243 10.829 11.305 11.797
beta_H[16,5] 9.861 0.182 9.492 9.862 10.212
beta_H[1,6] 10.118 0.232 9.749 10.089 10.678
beta_H[2,6] 11.470 0.110 11.253 11.470 11.684
beta_H[3,6] 10.617 0.158 10.283 10.624 10.905
beta_H[4,6] 12.900 0.876 11.165 12.909 14.569
beta_H[5,6] 5.652 0.618 4.553 5.614 6.934
beta_H[6,6] 8.745 0.565 7.483 8.781 9.770
beta_H[7,6] 10.544 0.584 9.404 10.550 11.703
beta_H[8,6] 9.518 0.288 8.911 9.549 9.963
beta_H[9,6] 8.311 0.777 6.819 8.306 9.875
beta_H[10,6] 9.137 0.340 8.399 9.156 9.733
beta_H[11,6] 10.676 0.369 9.895 10.711 11.316
beta_H[12,6] 9.378 0.242 8.931 9.370 9.882
beta_H[13,6] 11.089 0.152 10.802 11.086 11.399
beta_H[14,6] 9.920 0.285 9.350 9.926 10.463
beta_H[15,6] 10.535 0.435 9.676 10.542 11.366
beta_H[16,6] 10.576 0.254 10.044 10.585 11.061
beta_H[1,7] 10.152 1.230 7.058 10.431 11.873
beta_H[2,7] 12.107 0.519 10.950 12.150 13.055
beta_H[3,7] 10.309 0.650 8.926 10.350 11.489
beta_H[4,7] 1.547 4.531 -6.973 1.465 10.913
beta_H[5,7] 5.835 2.010 1.910 5.882 9.728
beta_H[6,7] 7.151 3.176 -0.289 8.272 11.132
beta_H[7,7] 7.688 2.950 1.872 7.694 13.573
beta_H[8,7] 10.893 1.108 9.071 10.851 12.866
beta_H[9,7] 4.835 4.009 -3.183 4.957 12.714
beta_H[10,7] 10.034 1.603 7.337 9.897 13.604
beta_H[11,7] 11.412 1.799 8.231 11.187 15.472
beta_H[12,7] 9.776 0.967 7.793 9.854 11.344
beta_H[13,7] 11.717 0.701 10.123 11.783 12.899
beta_H[14,7] 10.608 0.919 8.755 10.634 12.249
beta_H[15,7] 14.032 2.268 9.506 13.989 18.549
beta_H[16,7] 12.905 1.367 10.522 12.765 15.874
beta0_H[1] 8.229 18.801 -30.333 8.612 45.364
beta0_H[2] 10.603 7.994 -5.651 10.560 26.712
beta0_H[3] 9.614 10.255 -11.618 9.696 29.970
beta0_H[4] 2.641 185.320 -391.042 3.894 385.816
beta0_H[5] 4.578 26.327 -49.486 4.371 56.832
beta0_H[6] 5.411 53.251 -111.840 7.480 119.066
beta0_H[7] 1.635 139.169 -284.590 1.813 285.176
beta0_H[8] 7.200 35.536 -21.665 6.402 34.899
beta0_H[9] 5.968 125.030 -257.402 7.654 253.413
beta0_H[10] 7.595 38.175 -70.781 8.164 87.741
beta0_H[11] 10.020 51.905 -100.290 10.574 119.946
beta0_H[12] 7.089 12.760 -18.865 6.872 33.800
beta0_H[13] 9.847 11.038 -11.971 9.792 30.137
beta0_H[14] 7.073 12.204 -17.320 7.193 31.032
beta0_H[15] 12.335 115.618 -210.865 11.936 260.816
beta0_H[16] 7.599 32.332 -59.855 7.933 73.886